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KJNttECflN 2V!ftmATIQN SHBB0OT1HBS 30R REAL-TIME! PAHAIXSL QPERAriOK 

Consider a digital computer which operated on two accumulators &a 
parallel *» i.e., stores and manipulates a pair of words (a,bh 

la such real-time applications as conversion from polar coordi- 
nates {r,®) to cortesian coordinates (x,y} 9 or vie© versa,, time for 
calculation can now be halved by evaluating a different fraction in 
each side of the accumulator. 
For example; 

(i) (r t e)~*(x»y) vim 

Calculate jB * 2/tx©, given $ in radians t or^f * 0/<?o » 

given $ in degrees, put &.* fractional part of^$ and 

calculate !-«(*• 

Then via subroutine: (sin 1»«U sin * l ~) 

multiply by r (x,y) 

{it) {x 9 y)-**(r,8} vias 

if x » y, then r •» ^TS'X, <» 1/3*. 
if x«^y, calculate u<l where u «■ x/y (or y/x) 
Often via subroutine {l/B/fH^j{^f tL L}ru. ) 
Multiply by (y,u) : <r/fe,«0 

This last subroutine could be tab3je look-up with quartle inter- 
polation {see Method #lc) P or via summation of a truncated expansion 
for the function in a Jcurier-Ghebyshev series, using T {x) over 
{-1,1) or T *(x) over [0,1) {see Method #3)* 



DGL-1XS 

Meth od Jjl 

The r:jost straightforward laethoci is simultaneous table look-up 
and interpolation (car* Todd, Trans, Sy*ja>* App» Mh 4 , v2 f pp, 102-4)., 

( a ) Llicaa r Int e rpolation* 

iiiiii mu yn \ :i m ,tmit a t uM \ *n Mwn im»i— ■>— «» 

Store ^ - f(xj, i » 0»l,3*« tt and 5^* 'I&e calculation cycle 
is of ainliaom length, but to g\>t accuracy "to five significant figures » 
many tabular values must bo stored « Nearly all the program is a 
lengthy table , 30 the program 1 *ices too long to read in ana occupies 
too much apace in core joemory* 

(b) Quadratic interpolation* 

■Jtoi— * I ■ m i i ■ M m I iw wwy » i »i '■■ i i»i n i 

Store f< and J f* ami u& J- Everett's or Besael's interpolation 
formula* Fewer table values .'jeed be stored (only a smll percentage 
of those needed for (a) and fchis saving in length far outweighs the 
disadvantage of a longer calculation cycle,, 
(c) 

Store t^ and S m f t * S t ± ~0„ 18393 6V. This yields nearly 
fourth-order accuracy, so still fewer tabular values need be stored 
(leas than half the nuiaibai? for (h)) s yot the calculation cycle is 
exactly the same as in (*>}• (The only penalty is more preliminary 

f. is obtained only after 
calculating a f « and S f « ) 

Everett's formula! f p i= <^f a t ?$, + £«&%+£ ffii^Q$ U ' ° *' 
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u n-'/^-v\ a -^-tH /-*>»£ a * ills!) } £ 7 z/*Ji^: 

where p r 'M a ~~)> ^ * ( ^< * ' if r) ^o <g ' *~ ' £- 

Here o> ~ S' X J$$ } %-$(p)~- Ar ^~j£rH approximately constant over 
di f>^ I o (varies frc;a ~.15 to -»S0} o 
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